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Control of an atom laser using feedback 
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A generalised method of using feedback to control Bose-Einstein condensates is introduced. The 
condensates are modelled by the Gross-Pitaevskii equation, so only semiclassical fluctations can 
be suppressed, and back-action from the measurement is ignored. We show that for any available 
control, a feedback scheme can be found to reduce the energy while the appropriate moment is 
still dynamic. We demonstrate these schemes by considering a condensate trapped in a harmonic 
potential that can be modulated in strength and position. The formalism of our feedback scheme 
also allows the inclusion of certain types of non-linear controls. If the non-linear interaction between 
the atoms can be controlled via a Feshbach resonance, we show that the feedback process can operate 
with a much higher efficiency. 

PACS numbers: 03.75.Pp, 03.65.Sq, 05.45.-a 



I. INTRODUCTION 

In recent years, we have seen the first examples of the 
atom laser, a device similar to the optical lasei\providing 
a coherent, Bose-condensed output beam Q, IS 0> HJ • 
The development of the atom laser past the demonstra- 
tion stage, particularly the development of the pumped 
atom laser, is an important goal in atom optics. In many 
applications, it is the high spectral flux and coherence 
provided by a pumped laser that is critical. Pumping, 
however, is difficult to implement with atoms and can 
lead to classical noise that far exceeds the suppression 
of quantum noise, or line narrowing, that would be ex- 
pected from a pumped system. This paper presents a 
method of suppressing the classical noise on a pumped 
atom laser beam by feedback to the condensate, with the 
aim of achieving quantum noise limited operation. 

As with light, the matter waves from an atom laser 
can be coherently reflected, focused, beam split, and po- 
larised 0. These are the basic operations performed in 
all optics experiments and through these operations every 
linear, non-linear and quantum optics experiment has its 
analogue when performed with atoms. Although bosonic 
atoms and photons both exhibit Bose-stimulated scat- 
tering that is fundamental to laser operation 0, @ , there 
are significant and interesting differences. The free space 
dispersion relation for atoms leads to spatial broaden- 
ing of pulses in vacuum. Atoms interact with each other 
and display non-linear effects in the absence of another 
medium. Atoms display far more complex polarisation 
states, move slowly and can be readily produced with 
wavelengths much shorter than is available from an op- 
tical laser. These are ideal properties in many precision 
measurement and quantum information applications. 

The present state of the art in atom lasers is an un- 
pumped Bose-Einstein condensate (BEC) that serves as 
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a source for a propagating matter wave beam. Atoms 
are outcoupled from the condensate via an RF, or a 
Raman transition that coherently flips a trapped spin 
state to an untrapped state. There have been several ex- 
perimental investigations of the properties of atom laser 
beams. Both temporal and spatial coherence have been 
measured, and it has been demonstrated that RF outcou- 
pling preserves the coherence of the condensate 0, ^| . 
The beam divergence has been measured [l it . and there 
has been one real time measurement of the flux of an 
atom laser beam |l2j. The four wave mixing experi- 
ments performed by the NIST group were the first exper- 
iments to exploit the inherent nonlincarity of atoms in a 
controlled fashion and, furthermore, demonstrated that 
the Raman outcoupling process also preserved the coher- 
ence displayed by the condensate [t|. There have been 
two early experiments reporting squeezing in atom laser 
beams |l4l Il5j . Despite these pioneering experiments, 
there is a significant amount of development needed if 
the atom laser is to become a generally applicable and 
useful tool in quantum atom optics. 

High spectral flux in optical lasers is generated through 
a competition between a depletable pumping mechanism 
that operates at the same time as the damping. The 
linewidth of a pumped laser is much narrower than the 
linewidth of the cavity determined by the cavity lifetime. 
In a pumped laser, there is Bose enhancement of the 
scattering rate into the lasing mode resulting in line nar- 
rowing 
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The line narrowing, or suppres- 
sion of quantum noise associated with pumping an atom 
laser, is a very desirable but as yet unrealised property. 
Quantum field theory is required to calculate the quan- 
tum noise limited linewidth of an atom laser with inter- 
actions. Wiseman and Thomsen have studied the quan- 
tum noise on an atom laser beam outcoupled from a sin- 
gle mode condensate and have included feedback in their 
model. Atomic collisions turn number fluctuations into 
phase fluctuations significantly increasing the linewidth. 
A continuous QND feedback scheme can be used to can- 
cel this linewidth broadening [2(j. It would be difficult 
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to treat both quantum and classical noise in the same 
model, as the full quantum field theory is only tractable 
in the limit of a few modes, whereas the classical noise 
is intrinsically a multimode effect. There is no guaran- 
tee that a real atom laser would operate at the quantum 
noise limit, and it is likely that we must design pumping 
schemes very carefully and use feedback to approach the 
quantum noise limit. It is this goal that motivates the 
present work. 

The classical noise on a pumped atom laser can be 
studied with multimode semiclassical Gross-Pitaevski 
(GP) models j2j]]. In a recent paper, it was shown that 
an atom laser pumped by a non mode-selective pumping 
scheme was unstable below a critical value of the scatter- 
ing length leading to significant classical noise on the out- 
coupled beam [22]. It would seem sensible to adjust the 
scattering length via a Feschbach resonance to a suitably 
large value to stabilise the atom laser and reduce classical 
noise. Quantum and classical noise scale oppositely with 
scattering length, quantum noise increasing with scatter- 
ing length and classical noise decreasing. The solution is 
to operate at low scattering length and either use mode 
selective pumping to stabilise the laser, or to suppress 
classical noise by feeding back to the condensate. Mode 
selective pumping would appear to be difficult to imple- 
ment, and it is the second option, suppression of classical 
noise by feedback, that we investigate here. 

Any realistic feedback scheme will require a detector 
to measure classical noise, and a control to feedback to 
in order to suppress the motion of the condensate. The 
entire feedback loop must have sufficient bandwidth and 
must be minimally destructive. The design of minimally 
destructive detectors for real time measurement and feed- 
back to stabilise an atom laser was discussed in two recent 
papers 0, 113 ■ I n the present work, we have chosen to 
feedback to realistic controls provided by the magnetic 
trap to ensure straightforward implementation in an ex- 
periment. 



II. CONTROL OF A CONDENSATE 

The choice of an effective feedback scheme is largely 
determined by the available methods of controlling that 
system. For a BEC, these controls can correspond either 
to perturbations in the trap potential, or changes in the 
interactions between the atoms. We examine the feed- 
back scheme required to control a BEC in three dimen- 
sions in an arbitary potential. We model the system by 
the Gross-Pitaevskii equation. We assume that it is pos- 
sible to control a set of external potentials J^i a i(t)fi( r ) 
and spatially dependent nonlinear interaction strengths 
2j bj(t)gj(r) with time dependent amplitudes. With the 
feedback switched on, the equation of motion is: 

m^A^H^t) (i) 



with 

H = H Q + £o 4 (t)/ i (r) + ]T ^(^(r)!^ 2 (2) 

i j 

H = f+V (r) + UoH\ 2 , T=-^V 2 (3) 

Zm 

The a,i(t)'s and &i(i)'s are the set of controls used to ma- 
nipulate the potentials. We consider a condensate ini- 
tially evolving under the Hamiltonian Ho- Unless the 
system is initially in the ground state, we want to reduce 
the energy, given by: 

Eo^) = (f + Vo)+ l -(Uo\i>\ 2 ) (4) 

Where the angle brackets denote (q) = J 4>*qipd 3 r, and 
the integral is over all space. By switching the feedback 
on, and then switching it back off again at some later 
time, we will typically have altered the value of Eq. It 
is important to note that in the presence of feedback, 
Eq does not represent the instantaneous energy, but the 
energy that the system would have if the feedback were 
to be suddenly switched off at that time. The rate of 
change of Eq while the feedback is switched on is: 

Eq = Jrif+Vo^J + riT+Vo^J + ^-^m^r (5) 

Where the dot ' denotes differentiation with respect to t. 
Using equation in equation and the fact that H 
is Hermitian gives: 

Using the divergence theorem gives 

= /E«,/ t (r)(fVV-^VV)d 3 r 

i 

- ^ f^6 j5j (r)|^(fVV-^V 2 f)d 3 r 

j 

= -X)°*(*)</i( r ))-5E 6 iW<fi( r jW> ( 6 ) 

i j 

It can be seen that setting cij(t) = Cj(/j(r)) and bi(t) = 
Ui(gi(r)\ip\ 2 ) , where q and Uj are positive constants, so 
that 

2 • 2 

E = -Y,Ci{fi(r)) 4E%(fli( r )W 2 ). (7) 

i j 

will always reduce Eq while there are oscillations present 
in the appropriate moments of the condensate. This is 
an important result as it illustrates a general scheme to 
reduce the energy from the condensate depending on the 
available controls. In the following sections, we demon- 
strate applying this feedback scheme to particular ex- 
amples. In section IIIII we investigate how we can use 
feedback to control a linear (f/ = 0) system in a har- 
monic potential. In section llVl we demonstrate control 
of a Bose-Einstein condensate in a harmonic potential. 
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III. USING THE FEEDBACK SCHEME TO 
CONTROL A LINEAR HARMONIC 
OSCILLATOR 

A. Harmonic oscillator with linear controls 

We now consider the specific example of the linear 
(Uq = 0) Schrodinger equation in one dimension with a 
harmonic potential, ie: Vq = \x 2 (in harmonic oscillator 



units where x is measured in units of the length 




time t is measured in units of the time a-> _1 and energy is 
measured in units of Hu>, where to is the harmonic trap- 
ping frequency). We use as our controls the position of 
the minimum of the potential, and the strength of the 
potential. This system is a good model of a BEC in 
either a magnetic or an optical trap, which are both ap- 
proximately harmonic near the potential minimum, and 
can be modulated in intensity. The equation of motion 
is then 

iip= (f + V + a!{t)x + a 2 (t)x 2 )i} (8) 

Setting a±(t) = c\{x) and ai (t) = C2(x 2 } in accordance 
with our feedback scheme gives 

E = - Cl {x) 2 - c 2 (x 2 ) 2 (9) 

This will guarantee that Eq is always reduced while there 
are fluctuations in (x) and (x 2 ), but the rate can be opti- 
mized by carefully selecting the value ci and C2- We can 
calculate a dynamical equation for (x) using Ehrenfcst's 
theorem |2(j 

<i> = -(^r^) = -(! + 2« 2 (i))(z> - cxix) (10) 

This is mathematically identical to a classical damped 
harmonic oscillator. Critical damping will occur when 
ci = 2^/1 + 2<i2. The dynamic equation for (x 2 ) isn't a 
simple linear harmonic oscillator, so we found an appro- 
priate magnitude of C2 numerically. 

Equation [3] was integrated numerically using a pseu- 
dospectral method with a fourth-order Runge-Kutta time 
step using MATLAB. The feedback initially turned 
off, and then switched on at time t = 20. Figure shows 
how the oscillations in (x 2 ) are damped for different val- 
ues of C2- It appears that critical damping occurs when 
c 2 ~ 1, and this value will be used for all subsequent 
calculations. 

We next demonstrate that the two moments of feed- 
back can be used together to reduce energy from the 
system. Figure (j2J shows the system initially in a non- 
stationary state. The feedback is turned on at time 
t = 20, and oscillations in (x) and (x 2 ) are quickly re- 
duced. Eq is reduced until it is ^, which is the energy of 
the ground state wave function in a harmonic potential. 

In this particular example, the energy is reduced until 
the system is in the ground state. Equation [S] shows that 




t (Harmonic oscillator units) 

FIG. 1: Oscillations in condensate width versus time for (a) 
C2 = 0.05; (b) C2 = 1; (c) C2 = 5. It can be seen that (a) 
is underdamped, (b) is close to critical damping, and (c) is 
overdamped. (x 2 ) and t are measured in harmonic oscillator 
units. 
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FIG. 2: All quantities measured in harmonic oscillator units. 
Both modes of feedback working simultaneously on a system. 
The density profile of the initial condition is shown on the 
right with the solid black line, in comparison to the ground 
state density profile, indicated by the dashed line. The central 
density is the density at the point x — 0. The energy is 
reduced to Eo = 0.5, which is the ground state energy of the 
harmonic oscillator. 



the energy will only be reduced when there are oscilla- 
tions in (a;) and (a; 2 ), so once the system is in a state 
where (x) = and (a; 2 ) = 0, the feedback will no longer 
reduce the energy. Obviously, energy eigenstates will dis- 
play no error signal, but these are not a problem as they 
are single mode and all expectation values of observables 
display no time dependence. Using the harmonic oscilla- 
tor ladder operators (a = -?(^ +x), = — x)) 



4 



Density 



0.5 



0.4 



0.3 



0.2 



0.1 




1 

0.5 



0.01 


-0.01 

1.51 
1.5 
1.49 

1.6 
1.5 
1.4 



Central density 



x 2 > 



5 

x (H.O. Units) 



10 



20 



30 



40 



50 



t (H.O. Units) 



FIG. 3: A state with no oscillations in (a;) and (x 2 ). The 
feedback does nothing to reduce the energy as there is no 
error signal. 



we can write x = — at). In the absence of error 

signals {ai(t) = a,2(t) = 0), we can use the Heisenberg 
equation of motion to calculate (x) and (x 2 ). By setting 
these equal to zero, we get a condition for our zero error 
signal states 



n=0 



Vn + T(a* n+1 a n e lt + a* n a n+1 e lt ) = 



n=0 



Vn+ l\/n + 2(a 



n+2 UL n 



a n e 2it -a* n a n+2 e 2lt ) 



(11) 



(12) 



where |n) are the energy eigenstates (Ho\n) = E n \n}), 
and a n are their coefficients \ip) = Yl^Lo c^ne'^^^ln). 
This shows us that there are an infinite number of non- 
stationary states that display no error signal. 

This result demonstrates that feedback using these 
controls will not always be effective, as the system may 
be attracted to one of these states rather than an eigen- 
state. In these non-stationary states with no error signal, 
the energy will not be further reduced, and semiclassical 
fluctuations will continue. Figure© shows an example 
of such a state. It displays no oscillations in (x) and (x 2 ), 
and the feedback does nothing to reduce the energy. The 
oscillations in the density at the centre of the trap are 
included to demonstrate that the condensate is dynamic. 
Obviously, our two error signals are insufficient to reduce 
dynamics fluctuations for the system in general. Our 
choice of error signal is governed by the controls we have 
available to us. We chose the curvature and position of 
the minimum of the harmonic potential as our controls 
as they are easy to manipulate in current experimental 
traps. In the next section we introduce a time dependent 
nonlinear interaction in an attempt to produce a feedback 
scheme that will remove all the semiclassical fluctuations. 
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FIG. 4: A condensate in the same initial state as[3] but feed- 
ing back using a time dependent nonlinear interaction with 
iii = 5 as well as the two trap parameters. In this case the ad- 
ditional error signal allows the feedback to reduce the energy 
until it is the ground state energy. The condensate number 
was normalized to unity for this example. 



B. Harmonic oscillator with a nonlinear control 

It is possible to tune the non-linear interaction between 
atoms in a Bose-Einstein condensate by controlling the 
magnetic field close to a Feshbach resonance H|. In ex- 
perimental systems, this is equivalent to controlling the 
bias magnetic field in a magnetic trap, or applying a 
constant magnetic field in an optical trap, and it has 
been achieved with considerable finesse in many recent 
experiments 29] . Adding a time dependent interaction 
between the atoms gives the equation of motion: 

iip=(f + V + ai{t)x + a 2 (t)x 2 + bi (*)|V>| 2 )V> (13) 

Setting bi(t) = ui(\ip\ 2 ) in accordance with our feedback 
scheme will always reduce Eq. Figure^shows a system in 
the same initial state as figure [21 but with the additional 
control. The additional error signal allows us to perturb 
the system from the stable state, and the energy is re- 
duced to the ground state energy. We have demonstrated 
how we can use feedback effectively to remove energy 
from nonstationary states in the linear regime (Uo = 0). 
In the follow section we look at the more physically real- 
istic example of a Bose-Einstein condensate with a strong 
nonlinear interaction. 



IV. CONTROLLING A BOSE-EINSTEIN 
CONDENSATE WITH FEEDBACK 

We use as our next example the more realistic system 
of a Bose-Einstein condensate with strong interatomic 
interactions in a harmonic trap. We begin by just using 
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FIG. 5: Feedback on a condensate with a large nonlinear in- 
teraction (UO = 100, condensate number normalized to unity) 
using x and x 2 as our controls for the time dependent poten- 
tial. The density profile of the initial state is shown on the left 
with a solid line, compared to the ground state with a dashed 
line. Oscillations in (x) and (x 2 ) are reduced and the energy 
is reduced to E ~ 8.51, which is the ground state energy, ci 
was chosen to be 0.05 for this example. 



the two trap controls as described in section llTTl to reduce 
the energy. Figure shows a condensate that is initially 
in an excited state, and the two modes of feedback reduce 
the energy until it is in the ground state. This is a special 
case, however, and figureElshows the feedback acting on a 
more general initial state. The energy is quickly removed 
from the two controlled moments, but there is still energy 
left in higher energy excitations. In contrast to the linear 
system, the motion in these higher moments is coupled 
into the controlled modes via the nonlinear interaction, 
and hence slowly reduced. This is an inefficient process 
that may be allieviated by including the time dependent 
interaction strength as a third control. Figure0compares 
the results of using all three feedback controls on a BEC 
with non-zero interaction with the effects of using only 
the linear controls. The use of the non-linear feedback 
dramatically accelerates the energy removal process after 
the rapid initial control due to the linear controls. 



V. CONCLUSION 

We have described a feedback scheme for reducing en- 
ergy from a BEC in an arbitrary potential with an ar- 
bitrary set of controls. This reduces the semiclassical 
fluctuations in the condensate, a process that will be es- 
sential for producing high quality atom lasers. In the case 
of a linear harmonic oscillator with a modulated trapping 
potential, we demonstrated that energy can only be ex- 
tracted from the moments in the motion corresponding 
to the moments present in the available controls. The 
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FIG. 6: Feedback on a condensate with a large nonlinear in- 
teraction (Uo = 100, condensate number normalized to unity) 
in a different initial state. The feedback quickly removes en- 
ergy from the two controlled modes, but energy in higher 
order excitations is more slowly reduced as it is coupled into 
the controlled modes via the nonlinear interaction. 
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FIG. 7: Comparison of energy reduction by feedback 
with and without the time dependent nonlinear interaction 
strength. The Solid line is Eo for[fj| and the dashed line is Eq 
with the time dependent nonlinear interaction included for 
ui = 1000. 



ability to modulate the nonlinear interaction between the 
atoms provides a feedback scheme that can control a far 
greater range of inital states. Formally, any eigenstate 
will be unaffected by the feedback scheme, but as our 
scheme can only remove energy from the system, a slight 
perturbation will usually result in the system coming to 
steady state in a lower energy eigenstate. 

In the case of a Bose-Einstein condensate with a large 
nonlinear interaction, there is already coupling between 
different modes of oscillations. This means that each 
mode of feedback can remove energy from more than one 
mode of oscillation. This indirect method of extracting 
energy from the higher modes is quite inefficient. Adding 
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a nonlinear control improves the efficiency of the feedback modes, 
because it directly removes energy from a larger range of 
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